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Abstract 

The optical response of a coupled nanowire dimer is studied using a fully quantum mechan- 
ical approach. The translational invariance of the system allows to apply the time-dependent 
density functional theory for the plasmonic dimer with the largest size considered so far in 
quantum calculations. Detailed comparisons with results from classical electromagnetic cal- 
culations based on local and non local hydrodynamic response, as well as with results of the 
recently developed quantum corrected model is performed. We show that electron tunneling 
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and dynamical screening are the major nonlocal quantum effects determining the plasmonic 
modes and field enhancement in the system. Account for the electron tunneling at small junc- 
tion sizes allows semi-quantitative description of quantum results within classical framework. 
We also discuss the shortcomings of classical treatments using non-local dielectric permittiv- 
ities based on hydrodynamic models. Finally, the implications of the actual position of the 
screening charge density for the plasmon ruler applications are demonstrated. 
KEYWORDS Nanoparticle dimer, quantum plasmonics, field enhancement, plasmon ruler 

In metal nanoparticles, collective excitation of the valence electrons induced by an incident 
electromagnetic field, the localized plasmon, leads to plethora of optical phenomena of significant 
current interest. For instance, strong plasmonic enhancement of the local fields-""" opens a route 
to numerous practical applications, such as surface enhanced Raman scattering (SERS),-"- optical 
nano-antennas— 1^ allowing, for example, the control of radiation from single quantum emitters,— 
or generation of extreme ultraviolet pulses by non-linear high harmonic generation.— In nanopar- 
ticle assemblies, the hybridization of plasmonic modes can serve for guiding of the propagating 
fields,—!^ as well as it offers a way for rational engineering of desired optical response and local 
field profile.- Not only light harvesting and sensing properties can be thus greatly improved, but 
also the specific geometry dependence of the optical response can be used as plasmon ruler to 
determine the arrangement and nanoscale distances within chemical or biological species. ■^^~— 

Significant advances in fabrication and manipulation techniques allow nowadays for precise 
control of the geometry of the structure.—"— In particular, for nanoparticle assemblies the size of 
the gap between the adjacent nanoparticles can be brought below one nanometer so that the electron 
tunneling across the junction becomes possible. Thus, the plasmonic devices enter the quantum 
regime which represents a significant challenge for available theoretical approaches. Indeed, most 
of the descriptions of the optical response are based on the solution of the classical Maxwell equa- 
tions where the nanoparticles are modeled with sharp surfaces and the quantum nature of electrons 
forming the screening charge is neglected. Such effects as experimentally observed appearance 
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of the charge transfer plasmon for conductively coupled particles prior to direct contact,— or 
optical rectification— 1^ can not be addressed within a classical framework. As shown with re- 
cent quantum calculations—"^"— the main ingredients missing in the classical model are the spill 
out of conduction electrons outside the nanoparticle surfaces, and the finite spatial profile of the 
plasmon-induced screening charge. The use of a nonlocal dielectric functions-^^"^ can account 
for the latter effect, where the variation of electron density is smooth rather than infinitely sharp, 
as assumed in the classical local approach. However, the spill out of electrons outside the nanopar- 
ticle surfaces and associated tunneling across the narrow interparticle junctions requires special 
treatment. 

In this respect, the data available on the coupled nanoparticle dimer is quite revealing. This is 
a prototypical system for plasmonic nanoparticle coupling, which has been extensively studied.-"^ 
Theoretical studies based on the solution of the classical Maxwell equations predict a discontin- 
uous transition from the capacitively to conductively coupled particles. For vanishing junction 
size, the fields in the junction diverge and the plasmon modes experience diverging red shift as a 
result of the interaction between high charge densities induced at the opposite sides of the junc- 
tion.---'^"^ The charge transfer plasmon appears abruptly after the conductive contact.—"^ The 
nonlocal calculations based on the hydrodynamic model-^i^^i^^^"^ have shown that because of the 
finite spatial profile of the plasmon-induced screening charge, the fields in the junction stay finite 
albeit large. The number of plasmon resonances and their frequency shift is much reduced com- 
pared to classical predictions. On the other hand, it follows from the quantum treatments— >^~— 
that for narrow junctions, electron tunneling can short circuit the junction and quench the plasmon- 
induced field enhancement. The nanoparticles appear conductively connected prior to direct con- 
tact, and the transition between the non-touching and conductive contact regimes is continuous. In 
particular, the charge transfer plasmon associated with interparticle charge transfer—"— progres- 
sively emerges in the optical response of the system, as has been fully confirmed in recent experi- 
ments.—'^ These quantum effects could be reproduced with Quantum Corrected Model (QCM)^- 
that treats the junction between the nanoparticles as an effective medium allowing to account for 
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the quantum effects within the classical local Maxwell theory.— 

Present situation is thus characterised by the experimental ability to create the structures with 
sub-nanometer characteristic scales and several sometimes controversial theoretical approaches to 
address the optical properties of such a nanostructures. While quantum calculations provide an 
a priori exact answer to the problem, to date the application of rigorous quantum mechanical ap- 
proaches to coupled plasmonic systems have been limited to rather small systems and are extremely 
heavy numerically. It is thus of considerable interest to have a bench mark quantum results on pos- 
sibly large realistic system of coupled plasmonic particles allowing detailed comparison between 
full quantum and macroscopic Maxwell theories. This would allow to assess of the role of quan- 
tum mechanical effects and of the possibility to account for these effects in classical theory. In this 
article which concerns a strongly coupled nanowire dimer, we offer such a comparison. We show 
that the quantum mechanical results for the optical response of the system can be quantitatively 
reproduced with simple models. Local classical description of the system allows to account for the 
interparticle tunneling through the narrow junction and for the actual location of the dynamically 
induced screening charges. 



Model and computational aspects 

The system considered in this work consists of two identical nanowires in vacuum. The nanowires 
are infinitely extended along z-axis as sketched in Fig. [B Each nanowire has a circular cross- 
section of diameter D. The nanowires are separated by a variable distance S and the incident 
light is polarized along the inter-particle axis x. Coupled nanowires and nanorods are often used 
in SERS, nanoantenna, and plasmon guiding applications. 3,4,10,1 1^20^,53-57 7]^^^^ jj^jg system is 
well characterized both from experimental and theoretical point of view. Full numerical studies 
based on the solution of the classical Maxwell equation have been performed,— ^^i^"^ as well as 
semi-analytical studies based on transformation optics.-^^aiii^ The system has also been investi- 
gated using non-local hydrodynamic description.—"^ The high symmetry of the system with its 
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Figure 1: Sketch of geometry of the nanowire dimer. Two identical cylindrical nanowires are 
infinite along z-axis and have diameter D of the circular cross-section in the {x,y)- plane. The 
nanowires are separated by the junction of width S. The incident radiation is linearly polarised 
with electric field along x-axis. 

translational invariance along the z-axis allows us to address on the full quantum level the case of 
the cylinders with D = 9.8 nm. To our knowledge this is the largest size of the plasmonic dimer 
described so far within the time dependent density functional theory (TDDFT) framework. 

For the nanowires we adopt the cylindrical jellium model (JM). Despite its simplicity, this 
model captures the collective plasmonic modes of the conduction electrons and has demonstrated 
its predictive power for quantum effects in nanoparticle dimers.-^^*^ While obviously not providing 
chemical accuracy, the JM is well suited for the description of the nonlocal effects due to conduc- 
tion electrons such as dynamical screening of external field and tunneling as we discuss below. 
Along with the possibility to treat relatively large system on a full quantum level, the JM model al- 
lows for direct comparison between quantum and classical electromagnetic theory results. Indeed, 
the physics underlying the Drude model for the metal response as well as the refinement of the 
Drude model introduced by the nonlocal hydrodynamic model corresponds best to free electron 
metals. For noble metals, such as silver and gold, the contribution of the localised d-electrons to 
the screening would obscure the comparison. 

Within the JM, the ionic cores of the nanowire atoms are represented with uniform back- 
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ground charge density no = r^) . The screening radius rg is set equal to 4 ao (Bohr radius 



ao=0.053 nm) corresponding to Na metal. Sodium is a prototype system for which the JM per- 
forms particularly well in description of the finite size non-local effects on optical properties.— 
It should be emphasized that the qualitative conclusions drawn in this work are robust and inde- 
pendent of the particular choice of density parameter. Each of nanowires is characterised by the 
number of electrons A^e per unit length so that from the charge neutrality the nanowire diameter 
is D = 2^/4N^r[j3. The circle of diameter D provides the position of the jellium edge separating 
uniform positive background from the vacuum. The jellium edge is located at half a lattice con- 
stant a (a = 4.23 A for Na) in front of the last atomic plane at the surface. We have performed 
calculations for A^e = 40 and A^e = 100 with D40 = 6.2 nm and Dioo = 9.8 nm respectively. The 
Fermi energy in both cases is at 2.9 eV below the vacuum level. 

The quantum calculations of the absorption cross-section are based on the Kohn-Sham (KS) 
scheme of the time-dependent density functional theory (TDDFT).— We use the adiabatic lo- 
cal density approximation (ALDA) with the exchange-correlation functional of Gunnarson and 
Lundqvist.— Retardation effects can be neglected due to the small size of the system. Extended 
description of present numerical implementation can be found in Ref.^^ Therefore we only dis- 
cuss here specific aspects linked with present work on interacting nanowires. In brief, the Kohn- 
Sham orbitals of the ground state of interacting dimer Yj{x,y) are obtained from these of non- 
interacting distant cylinders by adiabatically reducing the separation. The xi/j{x,y) are discretized 
on the equidistant mesh in cartesian coordinates, and the time dependent Kohn-Sham equations 



are then solved with initial conditions ^'j{x,yj = 0) = In Eq. m is the free electron 

mass, and A = d^/ dx^ + d'^ / dy^. The effective potential of the system Vgff depends on the electron 
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density n{x,y,t) given by: 



nix.yj) = 2 £ Xj 1' • (2) 

j(Eocc. 

The sum runs over occupied states, the factor 2 results from the spin degeneracy, and Xj is the 
number of electronic states associated with z-motion along the nanowire. 

Xj = ^^J2{Ef-Ej), (3) 

where Ef is the Fermi energy, and Ej is the energy of the \irj{x,y) orbital. 

The dipole absorption cross-section per unit length is calculated from the density dynamics 
induced by impulsive perturbation as: Gabsiio) = ^^/m{a((u)}, with a{(o) being the dipolarpo- 
larizability (per unit length) of the system, and c - speed of light in vacuum. In order to analyse the 
temporal evolution of the system corresponding to the excitation of the given plasmonic resonance 
we perform the TDDFT calculations for the nanowire dimer subjected to an incident x-polarized 
laserpulse. Electric field of the pulse is given by =£'oexp — {^^) cosQ.t. The frequency 
of the pulse D. is set resonant with studied plasmonic mode. The duration of the pulse T (typically 
50 fs) is large enough for the narrow spectral width. The amplitude Eq is any small number to 
guarantee the linear response regime.— The snapshots of the induced densities, currents and fields 
presented lower in this paper are extracted at r ~ r/2. 

The classical electromagnetic calculations of the absorption cross-section have been performed 
with Comsol Multiphysics package, (version 4.2a, www.comsol.com). Local, nonlocal, and quan- 
tum corrected model description of the dielectric properties of the system were used. It is note- 
worthy that because of the small size of the plasmonic dimer the absorption cross-section closely 
correspond to extinction cross-section in this case. 

Within the local classical approach the dielectric constant of the nanowires is described with 
Drude model, which is a good approximation for the free electron metal as sodium considered in 
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this work. 

£(cq) = 1- , ^ . (4) 

where (Op is the bulk plasma frequency given by (Op = ^jAKno/m = 5.89 eV, and 7 accounts for 
the damping. 

The quantum corrected model (QCM) allows to introduce the tunneling effect into the classical 
Maxwell equations. We use the local approach where the nanowires are described with Drude 
model Eq. IH and the tunneling across the vacuum gap between the nanowires is accounted for by 
filling the junction with an effective dielectric medium.— It is described with the Drude model, 
similar to Eq. |4l 

£,ff{y,S, C0) = \- / (5) 

ft>(£t) + «7eff(j,5)) 

The effective damping 7eff models transition from resistive (large 5) to conductive (small 5) char- 
acter of the junction. Thus, it acquires a dependence on the separation distance S and y-coordinate. 
For very large S the QCM is exactly equivalent to local classical approach. No tunneling is possible 
and the vacuum gap limit is retrieved with Eeffiy.S, £0) = 1. In this work we use exactly the same 
parametrization of ggff as in the earlier publication, where the QCM is introduced and described in 
great detail. — 

Within the non-local hydrodynamic description the transverse component of the permittivity 
tensor is given by Eq. dH), and the longitudinal component acquires the wave vector k dependence: 

col 

£l(6)) = 1- . ^212 - 

(0[(0 + iy) - P^k^ 



The /3 parameter is given by /3 = a/3/ 5vf .-'^'^'^"^ In the present case, the Fermi velocity of 



conduction electrons Vf = ^3n^no/m = 1.05 10^ m/s, resulting in /3 = 0.81 10^ m/s. Provided 
transverse and the longitudinal components of the permittivity tensor the absorption cross-section 
is calculated with numerical approach as recently implemented by G. Toscano and coworkers 
within the Comsol Multiphysics package.— 
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Figure 2: Absorption cross section a for the single jellium nanowire of the diameter D = 6.2 nm 
(left) and D = 9.8 nm (right). Results are shown as function of the frequency CO of the incident 
radiation. The incoming field is the .x;-polarised plane wave. The TDDFT calculations are com- 
pared with results of the classical electromagnetic calculations using local (Drude) and nonlocal 
hydrodynamic response. (See the legend for definition of diffemt symbols used in the figure). 



Results and discussion 



Individual Nanowire 

We first characterize the linear optical response of individual nanowire. In Fig. [2] we show the 
dipole absorption cross-section o calculated with full quantum TDDFT method and with classical 
local Drude and non-local hydrodynamic approaches. The classical local approach with the Drude 
model given by Eq. |4]for the permittivity of the nanowire leads to the dipole resonance located at 
surface plasmon frequency (Ogp = cOp/Vl = 4.16 eV irrespectively of the nanowire diameter. We 
find that the absolute value of o and the width of the resonance as calculated with TDDFT is best 
reproduced with damping parameter 7= 0.275 eV for the D = 6.2 nm nanowire, and 7= 0.247 eV 
for the D = 9.S nm nanowire. The smaller damping obtained for larger nanowire indicates smaller 
decay due to surface scattering. Indeed, the surface scattering contribution to the damping can 
be expressed as 75 = Avp/Rd,^^ where = D/2 is the cylinder radius, and A is the system- 
dependent parameter. The calculated change of the damping is consistent with A = 0.5 in accord 
with available experimental and theoretical results for alkali clusters.— Obviously, through the 
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width of the single wire resonance the quantum size effects would affect the optical response of 
the nanowire dimer.— 

The dipole resonance calculated with TDDFT is red shifted from the classical (Ogp value but 
approaches (Ogp for increasing size of the system. For the D = 6.2 nm nanowire, the resonance is 
at (0 = 4.027 eV, and it is at (O = 4.072 eV for the Z) = 9.8 nm nanowire. Thus, present TDDFT 
results for cylindrical nanowire show the same trends as studied in detail for spherical alkali metal 
clusters.-^'^--2^The non-local effect at the origin of the red shift of the plasmon resonance is linked 
with dynamical screening of the fields by conduction electrons. The dipole resonance frequency 
can be found from:— "— 

ft)/G)sp = 1 -Re[d{0},^)]/R,^ + 0{R-^). (7) 

The real part of Feibelman's parameter Re[d{o))] gives the position of the centroid of the induced 
charge density.—"— When measured from the jellium edge, Re[d{o)s,^)] is positive for alkali metals, 
i.e. the screening charge is shifted into the vacuum because of the spill out of the conduction 
electrons outside the metal. From the calculated resonance frequency we obtain 7?e[J((Usp)] = 1 A 
in agreement with earlier calculations and experiment. '^^P^ ZL Importantly, for the noble metals 
such as Au and Ag the final size effects lead to the blue shift of the dipole plasmon resonance.—"^ 
The difference between alkali and noble metals is the localised d-electron contribution to the total 
screening in the latter case.— >^""— When the d-electron contribution is accounted for, Re[d{G))\ 
tums negative indicating that screening charge is predominantly induced inside metal.— "^^^^^ 

For the small diameter nanowire D = 6.2 nm, finite size effects appear in the TDDFT calcu- 
lations not only in the form of red shifts and line width changes, but also as clearly observable 
structures in the frequency dependence of the absorption cross section. These structures arise 
from the strong coupling of collective plasmon- and single electron-hole excitations. For large 
diameter nanowires, finite size effects become smaller and the plasmon resonance is much better 
defined. The TDDFT results also show shallow resonance in the optical absorption cross section 
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at ft) = 4.6 eV associated with a multipole plasmon (MP).— — The resonant frequency obtained 
with TDDFT is in good agreement with experimental data.— The density oscillations of the MP 
proceed within the layer of the spilled out charge as obviously can not be captured within the local 
classical theory. 

Finally, the green dotted lines in Fig. |2] show the absorption cross-section calculated with the 
classical non-local hydrodynamic model (NLHD).^^ In contrast to the TDDFT result and exper- 
imental data on alkali clusters,—"^ the NLHD model predicts a blue shift of the dipole plasmon 
frequency from COgp value. The reason for this is that independent of the metal, the plasmon- 
induced charges in NLHD are localized within the layer of the thickness P / COp below the metal 
surface. Thus, the effective Re[d^uiD{(0)] is always negative leading to the blue shift of the 
localised plasmon— in contradiction with the spill out effect known for alkali metals. However, 
negative i?e[JNLHD(ft))] inherent to NLHD description of conduction electrons and so the blue 
shift of the localised plasmon qualitatively coincides with full quantum result determined by the 
rf-electron screening for the noble metals.—"— 

Coupled Nanowires 

In Fig. [3] we show waterfall plots of the dipole absorption cross-section of the D = 6.2 nm and 
D = 9.8 nm nanowire dimers calculated with TDDFT. The results are presented as function of the 
frequency of the incident radiation for different widths S of the junction. The calculations have 
been performed for both positive and negative S, where the latter means a geometrical overlap 
of the nanowires. The 5 = case corresponds to the kissing cylinders-^^^ with touching jellium 
edges where the distance between the topmost atomic planes equals to the interlayer spacing, i.e. 
a continuous solid is formed at the contact point. 

For large positive separations the non-local effects including electron tunneling across the junc- 
tion are small, and the TDDFT results agree with earlier classical calculations— '^'^>^>^ (See also 
Fig. [5]). At 5 = 26.5 A the absorption spectrum is dominated by two resonant structures. First is 
the bonding dipole plasmon (DP) indicated with blue dotted line in Figure. It is formed from 
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Figure 3: Waterfall plot of the dipole absorption cross-section for the nanowire dimer in vacuum. 
The dimer consists of two Na nanowires of diameter D = 6.2 nm (left) and D = 9.8 nm (right) 
separated by variable distance S. The incoming field is the x-polarised plane wave. The centers of 
the particles are at x = ±(D + 5)/2), and S is negative for overlapping cylinders. S = —D would 
correspond to the limit of a single cylinder. TDDFT results are given as function of the frequency 
CO of the incident radiation for different separations S. For clarity a vertical shift proportional to 
the separation distance is introduced for each absorption spectrum. The red curves are used each 
5 ao ~ 2.65 A of S'-change. These are labeled with corresponding 5-values each 10 ao ~ 5.3 A 
of 5-change. The plasmonic modes responsible for the peaks in the absorption cross-section are 
labelled. DP stands for the Bonding Dipole Plasmon, QP for Bonding Quadrupole Plasmon, and 
CI for the lowest (dipole) Charge Transfer Plasmon, C2 for the higher energy Charge Transfer 
Plasmon. Vertical black dotted line shows the degenerate plasmon peak located close to cOsp. On 
the right panel the blue (green) dotted curve at 5 = 26.5 A shows results of the classical Drude 
(NLHD) calculations with adjusted parameters. Further details are give in the main text of the 
paper. 
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hybridization of the dipole plasmon modes of the individual nanowires. The second resonance 
indicated with black dotted line is formed by bonding quadrupole (QP) plasmon degenerated with 
higher order modes. It is slightly red shifted with respect to the surface plasmon frequency (Ogp- 
At 5 — 7- oo, the DP merges into the degenerate mode and the spectrum evolves into that of the indi- 
vidual nanowire. As the junction width S decreases, the DP shifts to lower frequencies because of 
the attractive interaction between the charges of opposite sign across the junction. Along with red 
shift of the DP, the QP (orange dotted line in Figure) splits from the degenerate mode, and shows 
the red shift with decreasing S. 

Despite the overall similarity of the results obtained with smaller and larger diameter nanowires 
there are some notable differences primarily caused by finite size effect. The resonances are much 
better defined for the larger D = 9. 8 nm nanowire. In this case the width of the resonances is smaller 
and the structures due to the interaction between the plasmon and electron-hole pair excitations 
disappear. For D = 6.2 nm nanowire, this interaction can even split the single plasmon peak into 
several lines. The red shift of the DP and QP modes with decreasing S is more pronounced for 
the Z) = 9.8 nm nanowire dimer. This is because of the scale (S/D) invariance of the plasmonic 
properties of nanostructure in the quasistatic limit. 

For junction widths below ~ 7 A, electron tunneling across the junction becomes important. 
The results obtained here for the nanowire dimer have much in common with those reported in 
quantum studies of touching spherical nanospheres,— and in recent experiments. ^^'^^ The DP 
resonance progressively disappears and the charge transfer plasmon mode (CI) emerges prior to 
the direct contact between the nanoparticles. CI appears as a broad shallow low-frequency peak 
at positive 5 ~ 1.5 A, and evolves into a well-defined resonance at 5 < 0. Similarly, because of 
electron tunneling, the QP mode continuously evolves into a higher order charge transfer plasmon 
mode C2 before direct contact between the nanowires. Thus, already at positive S the nanowires 
appear conductively connected showing characteristic charge transfer plasmon modes.— i^"— For 
a dimer with well established conductive contact, the CI and C2 modes experience a blue shift as 
also reported in classical calculations. -^^^ 
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Figure 4: Panels (a)-(k) Detailed analysis of the plasmon dynamics in the coupled Na nanowires 
of diameter D = 9.8 nm. The incident x-polarised laser pulse is at resonance with lowest (dipole 
plasmon mode of the system). Panels (a)-(j) present the snapshots of the induced charge density 
An, current density Jx, and field Ex for different junction widths S as indicated to the left of each 
row. The induced currents and fields are measured along interparticle x-axis. Positive (negative) 
values correspond to red (blue) color code. The induced densities are shown at the instant of time 
corresponding to the maximum dipole moment of the dimer. The induced currents and fields are 
shown at the instants of time when the induced fields in the junction reach maximum. Panel (k): 
Conductivity analysis. The current Jx measured on x-axis in the middle of the junction is plotted 
as a function of the normalized electric field at the same position. Different colors correspond to 
different separations S as explained in the insert. 
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To obtain further insights into the effect of the tunneling through the junction we have calcu- 
lated the electron density dynamics in the plasmonic dimer subjected to an incident x-polarized 
laser pulse resonant with lower (DP or CI) plasmonic modes. The panels (a)-(j) of Fig. |4] show 
snapshots of the induced charge density An, and x-component of the current density Jx, and field 
Ex for different junction widths S. Results are given for the case of Z) = 9.8 nm nanowires. The in- 
duced densities are shown at the instant of time corresponding to the maximum dipole moment of 
the dimer. The currents and fields are shown at the instant of time corresponding to the maximum 
induced field in the junction. 

For large separation 5 = 7.95 A the maximum induced dipole corresponds to the in-phase 
dipole polarisation of each nanocylinder expected for the DP mode. High charge densities are in- 
duced at the surfaces facing the junction resulting in large electric field enhancement l^o / Ein \ ~ 80. 
Here, Ein is the amplitude of the incident field and Eq is the amplitude of the field measured at the 
.x;-axis in the middle of the junction. The structure of the induced charges and fields are similar to 
previous classical results for coupled cylinders-^^^>^>^ and also resembles much to the case of 
metal sphere dimer.— The probability of tunnelling between nanowires is negligible and no current 
flows across the junction. Note that the maximum dipole polarisation corresponds to the instant of 
time when the maximum charge separation has occurred and the currents inside nanoparticles are 
minimal. For the reduced width S = 2.65 A, very similar profiles for induced densities and fields 
are obtained. Therefore, we do not show Ex but focus on the induced current. The junction width is 
now sufficiently small to allow weak electron tunnelling between the nanoparticles. Large optical 
field in the junction acts as large BIAS in the scanning tunneling microscope causing tunneling cur- 
rent across the junction.— Thus, the junction shows the resistive character with maximum current 
between nanowires reached at the maximum field and consequently at maximum induced dipole. 

Further reduction of the separation S increases the tunneling probability and short circuits the 
junction. When conductive contact is formed, the DP mode disappears and the CI mode emerges 
in the absorption spectrum. Panels (e) and (h) Fig. |4] show the induced charges for 5 = and 
S = —2.65 A, respectively. The results are very similar, i.e. already for the 5 = case of kissing 
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cylinders the charge transfer plasmon mode is well developed (see also Fig. |3l and Fig. |5]). The 
maximum dipole moment of the dimer corresponds to the oppositely charged nanowires. Max- 
imum currents and fields [panels (f), (g), (i), (j)] are reached when the total dipole moment of 
the system is minimum (compare with S = 2.65 A case). Precisely, this is large current flowing 
trough the entire system that builds the dipole polarisation with opposite charges at left and right 
nanowires consistent with resonant excitation of the CI mode. The field enhancement is about 30 
for both 5 = and S = —2.65 A separations. However, the fields are screened at the center of the 
junction, and the maximum fields are located at its sides. This is similar to the classical result for 
the overlapping cylinders.— 

Panel (k) of Fig. |4] provides further analysis for the evolution from resistive to conductive 
character of the junction with decreasing S. The current J^{t) on the x-axis at the middle of the 
junction between the nanowires is shown as a function of the field ^^(O at the same location. At 
large separations, the linear relation Jx{t) = gEx{t) between the current and the local field shows 
that the junction is resistive. The increase of the slope g when S is reduced from 5 = 5.3 A to 
5 = 2.65 A is because of the increase of the tunneling probability. For 5=1 .06 A and 5 = A the 
Ex{t) and Jx{t) acquire a relative phase. Since the field envelope of the incident pulse grows in time 
for the time interval shown in Fig. IHk), the clock-wise rotation of the /c(?)[£',^(?)] curve implies 
that the current is retarded with respect to the field. The junction becomes conductive which is 
particularly apparent for 5 = 0. 

The TDDFT results show that for small junction widths the optical response is determined 
by electron tunneling. For large width of the junction the tunneling is negligible. Each cylinder 
responds on the self consistent field created by the incident radiation and the neighbour. The 
effects of nonlocal screening discussed for the single cylinder should a priori influence the optical 
response of the dimer as suggested with NLHD calculations—"^ and as inherently build in the 
TDDFT results. While the effect of tunneling is clearly seen with TDDFT calculations, revealing 
the role of the nonlocal screening requires comparison with classical results, in particular for the 
energies of the plasmonic modes. In this respect, present system offers a unique opportunity to test 
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Figure 5: Comparison of the optical properties of a nanowire dimer obtained with the full TDDFT 
calculations, with the quantum corrected model (QCM), with a classical Drude electromagnetic 
calculations (Classical), and with calculations based on the nonlocal hydrodynamic model (Non- 
local). The dimer consists of two Z) = 9.8 nm Na nanowires in vacuum. The incoming plane wave 
is polarized along the dimer axis x. Upper panels: Waterfall plots of the dipole absorption cross- 
section as function of the width of the junction S. Red curves correspond to separation distances 
of 5 = -5.3 A, -2.65 A, A, 2.65 A, 5.3 A, 7.95 A, and 10.6 A. For further details see caption 
of Fig 13] Lower panels: Color plots of the local field enhancement at the center of the junction for 
positive separations. Results are shown as function of the frequency CO of the incident radiation 
and separation distance S. The color code is explained at the bottom of the corresponding panels. 
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the validity of different classical models against full quantum results for relatively large system, 
with fully developed plasmonic modes, and with dielectric properties well described within free 
electron approximation. The rest of the paper is devoted to quantum vs classical comparison with 
particular focus on possibility to account for quantum effects within classical approach. 

Quantum vs classcal approaches 

Fig.[5]shows the present TDDFT results, the QCM results, results from classical calculations based 
on the Drude description of Na dielectric function, and results of the NLHD calculations. Waterfall 
plots of the absorption cross-section (upper panels) and contour plots of the field enhancement on 
the X-axis at the middle of the junction are presented as function of the frequency of the incoming 
plane wave and of the junction width S. Since for the individual nanowire TDDFT gives a red 
shift of the dipole plasmon resonance from the classical (Ogp position, while NLHD leads to a blue 
shift (Fig. |2l), we adjust the parameters of different models. For the sake of comparison of the 
junction width dependence we enforce an agreement between the TDDFT, QCM, classical Drude, 
and NLHD at the largest S calculated here with TDDFT, S = 26.5 A (see Fig. [3]). For the D = 9.8 
nm nanowire dimer we have used plasma frequency (Op = 5.8 eV and damping 7 = 0.247 eV 
in the QCM and classical Drude calculations. Thus, only small correction to the bulk plasma 
frequency (nominal value 5.89 eV) is required in this case. The NLHD model needs for much 
larger correction to compensate for the blue shift inherent to this treatment. We used (Op = 5.5 eV 
and 7= 0.16 eV. Adjustment of the nonlocality parameter /3 is also possible, but then the model 
looses its predictive power. Indeed, as far adjustment is limited to 0)p and 7 parameters, the finite 
size effects and so the need for correction will disappear with increasing size of the nanoobjects. 

The main features of the quantum results have been discussed in connection with Fig.[3]so here 
we will focus on comparisons between different model approaches. As follows from Fig. |5l the 
QCM does an excellent job in describing the TDDFT results over the entire range of separations 
S addressed here. The important features such as: the number of resonances; their distance de- 
pendence; the transition from the separated to conductively coupled regime are well reproduced. 
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In particular, in sheer contrast with classical theory— '^^'^ij^'a^ the change of the spectrum at the 
moment of contact is progressive and fields at the middle of the junction are quenched, not diverg- 
ing with decreasing junction size. Besides qualitative aspects, we find that the TDDFT and QCM 
agree semi-quantitatively as is further stressed in Fig. |6] This figure zooms into the most delicate 
interaction regime corresponding to the transition from the separate to overlapping nanowires. 

By construction of QCM, it is equivalent to the classical local Drude description for large 
positive S where tunneling is negligible. The good agreement with TDDFT data suggests that the 
pure local classical description is reasonable for a large S. At the same time the classical description 
fails at small 5 ~ 5 A, i.e. typically at two lattice constants between surface atomic planes that 
define the junction.— The accumulation of charges on the opposite sides of the junction leads 
to exaggerated coupling between nanowires resulting in diverging fields and too large number of 
resonances. Similarly, for negative S, the sharp edges of the junction, which are otherwise smeared 
out by the electron tunneling also result in too many hybridized resonances.— 

The NLHD description is free from the S = divergence problem.-^^di. The number of reso- 
nances remains small and the fields in the middle of the junction stay finite albeit too large. This 
can be easily understood thanks to the elegant transformation optics approach developed for in- 
teracting cylinders by Femandez-Dommguez and coworkers.— At positive S an analytic solution 
is found which depends on the renormalised parameter S ~ S + 2d, where 8 ^ ^ / (Op (1 A in 
the present case). The physics behind this shift is the localisation of the plasmon-induced charge 
below the surface as inherent for the NLHD model. As we discussed for the single cylinder, 
in terms of the dynamical screening theory, 5 is equivalent to an effective Feibelman parameter 
d = —Re[d^uxD{0))\- Thus, even for 5 = 0, the induced charges at the opposite sides of the junc- 
tion are actually separated by a finite distance 25. However, since the tunneling is not accounted 
for, the NLHD description fails to reproduce the quenching of the field enhancement at the middle 
of the junction for small positive S. The NLHD also fails to smoothen out the transition from sep- 
arated to overlapping regimes and gives an abrupt nonphysical transition. The number of modes is 
smaller than in the classical Drude description, but still larger than what is obtained in TDDFT or 
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Figure 6: Detailed comparison between TDDFT and QCM calculations. The dipole absorption 
cross-section of the D = 9.8 nm Na nanowire dimer is shown for small separations S corresponding 
to the strong tunneling regime and transition from separated to conductively coupled nanowires. 
The frequency range is zoomed at the transition from the bonding dipole (DP) to the lowest charge 
transfer (CI) plasmon. Waterfall plots of the optical absorption cross-section are shown for the 
separation distances changing from S = —2.65 A (lowest blue line) to 5 = 4.77 A (upper black 
line) in steps of 1 ao (0.53 A). For further details see caption of Fig. |3] 
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QCM calculations. 

Our final remark concerns the QP plasmon that evolves into the C2 charge transfer mode for 
negative S. The associated resonances are much less pronounced in the TDDFT cross-section 
compared to model approaches. One possible reason is finite size effect, where the the system size 
is not large enough for the corresponding density oscillations to be completely formed. However, 
the similarity between the TDDFT results for D = 6.2 nm and D = 9.8 nm nanowire dimers as 
shown in Fig. |3] suggests that this quantum size effect is small. We thus tentatively attribute the 
weaker high order resonances as obtained in TDDFT to the effect of the smearing of the induced 
surface charge densities. This a priori reduces the coupling between the dipole and higher order 
modes and consequently the intensity of the QP resonance in the dipole absorption cross-section. 

Dynamic screening 

With results shown in the previous Subsection, tunneling determines the optical properties of the 
system at small junction widths. For large width S tunneling is absent and nonlocal dynamic 
screening influences the optical properties. Similar to the case of the individual nanowire, the 
major nonlocal effect is the red shift of the frequency of the bonding dipole plasmon as compared 
to the classical prediction. However, for the dimer, this effect is much stronger, and does not 
disappear with increasing radius of nanoparticle. 

Fig. m presents the analysis of the role played by dynamic screening in determining the fre- 
quency of the bonding dipole plasmon. Panel (a) gives schematic representation of the location 
of plasmon induced screening charges in the junction of width S. Within the local classical ap- 
proach the screening charges are localised at geometrical surfaces of the cylinders separated by S. 
We recall that we use the definition of the geometrical surface such that it coincides with jellium 
edge. Within TDDFT, the real part of the centroid of the induced charge density Re[d{co)] corre- 
sponds to the position of the plasmon-induced surface charges with respect to the jellium edge of 
each cylinder.-2i"i^i2Mi For the present case of sodium nanowires, Re[d{(o)] is positive ^0.9 A 
in the frequency range of interest. This means that the screening charge is shifted by 0.9 A into 
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Figure 7: Dynamic screening, (a) Schematic representation of the location of plasmon induced 
screening charges in the junction. Within the local classical approach the screening charges are at 
geometrical surfaces of the cylinders (here equivalent to the jellium edges) separated by the junc- 
tion of width S. Within the TDDFT, the centroids of the screening charges (red areas) are located 
at Re[d{(o)] in front of the jellium edges and separated by S — 2Re[d{(0)]. In NLHD approach 
the centroids of the screening charges (blue areas) are located at distance 5 below the geometrical 
surface and separated by S + 25. (b) Energy of the dipole plasmon resonance as function of the 
junction width S. Dots: TDDFT results obtained for nanowire dimers formed by D = 6.2 nm and 
D — 9.8 nm nanowires (see the legend). Solid and dashed lines show results of classical Drude 
calculations for D = 6.2 nm (blue) and D = 9.8 nm (red) dimers. Dashed lines: separation S is 
measured between the jellium edges. Solid lines: separation S is measured between the centroids 
of the induced charges. For more details see the main text, (c) Energy of the dipole plasmon 
resonance as function of the junction width S. Dots: TDDFT results obtained for D = 9.8 nm 
nanowire dimer. Solid and dashed red lines show results of classical Drude calculations [same as 
(b)]. Dashed and dotted gray lines show results of NLHD calculations. Dashed line: separation S is 
measured between the jellium edges. Dotted line: separation S is measured between the centroids 
of the induced charges, (d) Same as (b), but results are shown as function of the scaled separation 
S/Rch where R^ = D/2 is the nanowire radius. 
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the vacuum and located at ^ 3 A outside the surface atomic plane of the nanoparticle. The dis- 
tance between the plasmon induced charges across the junction is thus E = 5 — 2Re[d{(o)] . As to 
the NLHD approach, with present choice of the non-locality parameter /3 it places the screening 
charges at 5 = 1 A inside the geometrical surface. Thus, within the NLHD approach the screening 
charges are separated by the distance 5 + 25. 

This insight provides an explanation for the junction width dependence of the DP modes shown 
in Fig. Hb) for two different nanowire dimers. While the overall agreement between the classical 
Drude (dashed lines) and the TDDFT (dots) calculations is good, the TDDFT result is system- 
atically redshifted. The difference is increasing with decreasing S indicating that the classical 
calculations underestimate the coupling between the dipole plasmons localized on each nanowire. 
Our finding that the plasmon-induced charges in the TDDFT calculations are outside the nanopar- 
ticle surfaces explains this effect. Indeed, the actual energy of the DP is given by the interaction 
between the screening charges across the junction. Thus, the TDDFT data obtained for junction 
width S should be compared with classical results calculated for E = 5 — 2Re[d{(o)] . If the DP dis- 
persion is plotted as a function of effective junction width given by the distance between plasmon 
induced screening charges, the agreement between TDDFT and the classical simulations becomes 
excellent. We note in passing that the step-structure of the TDDFT data for the smaller D = 6.2 nm 
dimer steams from the strong coupling between plasmon and electron-hole pair excitations. 

A similar analysis can be performed for the NLHD model calculations shown in Fig. ITJ^c) for 
the case of D = 9.8 nm nanowire dimer. The NLHD results appear at even higher frequencies than 
the classical Drude calculations, in even stronger disagreement with the TDDFT than the classical 
model. The explanation for this is that NLHD artificially places the plasmon-induced screening 
charges at a position ~ 5 = 1 A inside the metal surface. Thus, the separation between screening 
charges 5 + 25 in the NLHD model is larger than separation 5 in pure local theory, and too large 
compared to the quantum result E = S — 2Re[d{(0)]. This decreases the interaction and results in 
blue shifted DP modes compared to the classical model and TDDFT. If the DP dispersion is plotted 
as a function of effective junction width, agreement between NLHD, classical Drude, and TDDFT 
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results is obtained. 

The above results show that the TDDFT dispersion of the DP can be fully retrieved with local 
classical calculations. The main issue here is the size of the junction for which the calculations 
have to be performed. For geometry of the nanostructure, it seems convenient to define the surfaces 
of the objects as given by the jellium edges. In this case the width 5 = of the junction would 
correspond to continuous solid formed at the contact point, which is physically sound. On the 
other hand, classical calculations have to be performed for the effective junction width given by 
the actual separation between the screening charges. These conclusions have direct implications 
for the ultimate limit of resolution of plasmonic rulers.— 

□ □ 

For the free-electron Na surface, the screening charges are located at approximately 0.9 A (3 A) 
outside the jellium edge (surface atomic layer) meaning that for a Na-Na junction, the effective 
junction width would be 1.8 A (6 A) smaller than the physical junction width measured between 
the jellium edges (surface atomic planes). This conclusion holds not only for the nanowire case 
considered here, but for any junctions between Na surfaces. For silver and gold, analysis of the data 
on the blue shift of the dipole plasmon resonance in small clusters— i^i^^*^^'^ with Eq.|7]places the 
effective screening charges inside the jellium edge at 1 .5 ^0.85 A for silver, and at approximately 
1 .5 A for gold. Therefore, for an Ag-Ag (Au-Au) junctions, the effective junction width would be 
by 1.7 ^ 3 A (~ 3 A) larger than the physical junction width measured between the jellium edges 
and close to the junction width masured between the surface atomic planes. 

The use of the plasmon ruler relies on the universal dependence of the DP frequency on the 
scaled separation.—*^ In Fig.|7td) we show the TDDFT and classical results for the DP frequency 
of the D = 6.2 nm and D = 9.8 nm nanowire dimers as function of the scaled separation S/Rc\. 
The TDDFT data for both nanowire dimer sizes nearly falls on the unique curve provided that the 
separations S are sufficiently large that no tunneling occurs. This holds for S larger than typically 2 
lattice constants between surface atomic planes, which sets lower bound for the distances that can 
be actually measured with plasmon ruler. 
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Summary and Conclusions 

In conclusion, we have presented fully quantum mechanical study of the optical response of the 
plasmonic dimer formed by realistic size cylindrical nanowires. This system is also representative 
of interacting nanorods and is of relevance for SERS, plasmon ruler, and plasmon transport appli- 
cations. Translational invariance allows to apply the time-dependent density functional theory for 
the plasmonic dimer of largest size considered so far in quantum calculations. 

We considered nanowires made of Sodium which is a prototype for free-electron metal so 
that jellium model applies. The free-electron character of Na valence electrons implies that its 
permittivity can be well described with Drude model, as well as it is consistent with approximations 
behind the hydrodynamic approach to model the nonlocal character of the dielectric function. 
Thus, this is the system of choice allowing to set the full quantum TDDFT benchmark results, and 
to use these results to test different theoretical approaches addressing plasmonic response of the 
strongly coupled objects. This was one of the central goals of the present work. 

We have found that for the small junction widths, the optical response is determined by the 
quantum tunneling of conduction electrons across the potential barrier separating the nanowires. 
A decreasing junction size leads to progressive attenuation of the plasmon modes of separated 
nanowires and the emergence of charge transfer plasmon modes of conductively coupled dimer. As 
this happens, the fields in the middle of the junction are screened. The maximal field enhancement 
moves from the middle to the external regions of the junction. In this distance range the classical 
local Drude and nonlocal hydrodynamic model descriptions fail since they do not account for 
tunnelling. At variance, the quantum corrected model reproduces the TDDFT findings on a semi- 
quantitative level. 

For large separations S, the tunneling is negligible and the overall agreement between TDDFT, 
classical and the QCM results is good. Thus the QCM performs well over the entire range of 
separations studied here. For large S the agreement between the classical and TDDFT results can 
be further improved by taking into account the shift of the plasmon-induced charge density with 
respect to the geometrical surface of the nanoparticle. The latter can be defined as the jellium edge 
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as in the present work, or as the top most atomic layer at the surface. The effective junction width 
is then given by the separation of plasmon indicted charges at the opposite sides of the junction. 
Indeed, this is the interaction between these charges that determines the hybridization and energies 
of the modes. Introducing a simple distance correction into the classical calculations allows a full 
account of this non-local effect. This result has implications for the plasmon ruler concept and 
shows that care should be taken with respect to the definition of the separation which is actually 
measured. 

We hope that the results presented here contribute to the understanding of the role of quantum 
nonlocal effects in strongly coupled plasmonic systems, and help in elaborating efficient theoretical 
approaches for their prediction. 
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